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Abstract. 

We present new calibrations of the widely used H„ , [Oil] , and UV luminosity vs. star formation rate (SFR) 
relations. Using our evolutionary synthesis code GALEV we compute the different calibrations for 5 metallicities, 
from 1/50 solar up to 2.5 solar We find significant changes in the calibrations for lower metallicities compared to 
the standard caUbrations using solar input physics. 
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. 1. Introduction 

■ Determining the Star Formation Rate (SFR) of galaxies is one 

■ of the most important steps to understand their nature and evo- 
lution. This is usually done by using the luminosity of the 
Ha or the [Oil] line, as well as the luminosity at 1500A and 
2800A in the UV, as tracer of the ongoing star formation (SF) 
(Kennicutt 1998, Gallagher et al. 1989, Madau et al. 1998, 

. hereafter K98, G+89, M+98). These methods use the fact that 

I the UV continuum and the output of ionising Lyman contin- 
uum (Lyc) photons, responsible for the gaseous emission, are 
dominated by young massive stars with masses > 10 and 

. lifetimes < 2 x 10^ yr. 

Conventionally, the standard calibrations are derived on lo- 
cal samples of normal, i.e. big galaxies, assuming solar metal- 

] licity, as appropriate. In the local universe, most of the actively 

■ star-forming galaxies are moderate to low luminosity late-type 
and dwarf galaxies having subsolar metallicities, and a sig- 
nificant amount is contributed by these systems to the local 
SFR density (e.g. Brinchmann & Ellis 2000, Brinchmann et 
al. 2004). Going to higher redshift, spectroscopic studies are 
clearly biased towards the brightest and most metal rich sys- 
tems at every redshift. Even those, however, reveal substan- 
tially subsolar metallicities (e.g. Mehlert et al. 2003) and multi- 
band photometry of deep fields reaching the bulk of the intrin- 
sically fainter galaxy population will clearly be dominated by 
subsolar metallicity galaxies and protogalaxies. 

While the eff'ects of the initial mass function (IMF) (K98) 
or the absorption by dust are often discussed as sources of un- 
certainty in determining the SFR (e.g. Inoue et al. 2001), metal- 
licity efi'ects are rarely addressed. Concentrating on the SFR of 
star-bursting dwarf galaxies derived from H^^ line fluxes, the 
metallicity dependence was examined e.g. by Weilbacher & 
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Fritze - v. Alvensleben (2001). For a first extensive observa- 
tionally based investigation into the metallicity dependence of 
[Oil] as a SFR indicator see Kewley et al. (2004). 

We here investigate the effect of metallicity on the calibra- 
tion of the various SFR estimators, using our evolutionary syn- 
thesis code GALEV, described in Bicker et al. (2004), extended 
to include the gaseous emission, for 5 metallicities from 1/50 
solar up to 2.5 solar. 

We first present some details of the gaseous emission as in- 
cluded into our code in Sect. |2] We then discuss the derived 
calibrations aacounting for the effects of metallicity as well as 
their dependence on the mass limits of the IMF and on the stel- 
lar evolutionary input physics in Sect.|3] Finally we summarize 
our results in Sect.|5 



2. GALEV: galaxy evolution models 

Our chemically consistent evolutionary synthesis code GALEV 
is based on a modified version of Tinsley's equations for 
the chemical enrichment of the inter stellar medium and on 
isochrones from the Padova group (Bertelli et al. 1994) in the 
Nov. 1999 version, and the spectral library from Lejeune et al. 
(1997, 1998) for the spectral and photometric evolution of the 
stellar component. For a detailed description see Bicker et al. 
(2004). Now we have included the effects of gaseous emission, 
in terms of lines as well as continuous emission. This was al- 
ready included into our GALEV models for single bursts sin- 
gle metallicity stellar populations (SSPs) by Anders & Fritze 
- V. Alvensleben (2003), so we follow the method used there 
to implement the gaseous emission into our chemically consis- 
tent galaxy models. To clarify the effects of metallicity we use 
galaxy models with fixed metalhcities in this Letter. 
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2.1. Gaseous emission 

The gaseous emission is dominated by very hot stars, which 
produce hydrogen ionising (Lyc) photons. This means, that the 
emission is important in the early phases of an SSP (cf. Anders 
& Fritze - v. Alvensleben 2003), or, in case of a galaxy, in 
stages of (high) star formation activity. For every star on every 
isochrone we calculate the flux of Lyc photons from up-to-date 
non-LTE expanding model atmospheres (Schaerer & de Koter 
1997, Vacca et al. 1996, Smith et al. 2002). Summing up aU Lyc 
photons of all stars on all isochrones at every timestep gives us 
the total number of Lyc photons (NLyc) at each time. On the 
basis of NLyc we calculate the gaseous continuum emission and 
the hydrogen line fluxes as described in Kriiger et al. (1995) 
and Weilbacher et al. (2000). So the flux is given by: 

F(H^) = 4.757 X 10"" • NLyc • f, 

where f describes the fraction of Lyc photons actually involved 
in ionising the ISM. For higher metallicities (Z > 0.008) we 
assume that 30% (f = 0.7) of the Lyc photons are absorbed 
by dust immediately and cannot ionise the gas, i.e. f - 0.7 
(Mezger 1978, Weilbacher et al. 2000). For the lower metal- 
Ucities we take f = 1 . Note that in this respect our models are 
not completely self-consistent. The input we use implies that 
the most massive stars at low metallicity rapidly emerge from 
their dust cocoons. Lines other than are calculated from line 
ratios relative to Hyj. Line ratios for hydrogen lines are taken 
from Stasinska's (1984) theoretical models, while for lines of 
other elements we prefer to use empirical line ratios as derived 
from observations by Izotov et al. (1994, 1997) and Izotov & 
Thuan (1998). We include a total of 31 hydrogen lines up to the 
Brackett series and 36 lines of other elements. For the continu- 
ous emission and a detailed description see Anders & Fritze - v. 
Alvensleben (2003). We recall that the metallicity dependence 
of the isochrones induces a metallicity dependence for the flux 
of ionising Lyc photons that, in turn, leads to a metallicity de- 
pendence of the hydrogen hues. For lines of other elements, 
additional metallicity dependences come into play in a variety 
of ways. Line ratios themselves depend on the metalhcity of 
the gas in a specific and complicated way for each element/line 
and difl'erences in the typical electron temperatures and den- 
sities in difl'erent metalhcity envirormients add to this. That is 
the reason why we use for these Unes ratios derived from ob- 
servations of a large number of difl'erent galaxies at each of our 
subsolar metalhcity intervals. 

In addition to young massive stars white dwarfs also have a 
significant output of Lyc photons. Their emission, however, is 
associated with the planetary nebula phase, which is fairly short 
in comparison with our minimum timestep (~ 25 x 10^ yrs vs. 
4 X 10*" yrs). Hence, we decide to ignore the contribution of 
the WDs to the total Lyc flux of actively star-forming galaxies. 
Binary stars are not yet included in our models at that stage. 

2.2. ZAMS extension of the Isochrones 

The youngest age provided by the Padova isochrones that we 
use is 4 Myr. This corresponds to an upper mass limit of 
~ 70 Mq. For the gaseous emission, however, the youngest and 



Table 1. Calibration constants for the and [Oil] luminosity 
vs. SFR relations. (Values in brackets are for f = 1) 



M 



up 



0.0004 
0.004 
0.008 
0.02 
0.05 



100 Mo 



C[oii] 



10« 



erg s 



Mq yi"^ 



2.7 
2.3 

1.4 (2.1) 



0.5 
1.4 

1.5 (2.2) 



1.3 (1.8) 1.3 (1.9) 
1.2 (1.7) 1.2 (1.8) 



120 Mp 



Ch. 



C[oii] 



Mq yr- 



3.0 
2.6 

1.7 (2.3) 



0.5 
1.6 

1.8 (2.6) 



1.6 (2.2) 1.6 (2.3) 
1.4 (1.9) 1.4 (2.0) 



most massive stars play a very important role, because the num- 
ber of Lyc photons (NLyc) increases dramatically with stellar 
temperature, NLyc ~ T^^j. The original isochrones hence un- 
avoidably implied a severe underestimate of the Lyc flux. 

Stellar evolutionary tracks from the Padova group, how- 
ever, start from zero age main sequence (ZAMS) and are avail- 
able for stars up to 120 M©. Hence we decided to supplement 
our set of isochrones by adding ZAMS isochrone. Stars on the 
ZAMS are not evolved, so there is no need to consider equiva- 
lent stellar evolutionary stages, usually a challenge in convert- 
ing tracks to isochrones. Now we are able to go to higher upper 
masses than the original isochrones provided. 

For our models we now adopt an upper mass limit of 100 
Mq as our standard and investigate 120 Mq for comparison. 

3. Results 

To study the metallicity dependence of the various SFR trac- 
ers we calculate models with constant SFR at the 5 available 
metalhcities. We assume a Salpeter IMF with lower and upper 
mass limits of 0.15 and 100 Mq as our standard and explore 0.1 
and 120 Mq for comparison. We have also explored difl^erent 
star formation histories (SFH) and found no difl'erences in the 
calibrations as long as the SFR is evolving smoothly. For a dis- 
cussion of rapidly changing SFRs, like in the short starbursts 
in dwarf galaxies, see Weilbacher & Fritze - v. Alvensleben 
(2001). 

The derived calibrations for SFRs in terms of emission lines 
directly depend on the fraction f of actually ionising Lyc pho- 
tons. For comparison purposes, we also give the values for 
f = 1 in the case of higher metallicities. The calibrations are not 
corrected to any amount of dust. Such a correction can easily 
be applied by calculating the ratio between the flux emerging 
from a dusty galaxy Fdust and the corresponding unextincted 
flux Fq. For any extinction law k(A) this is given by 



Fdust (/I) 
Fo(^) 



= 10 



-0.4-E(B-V) 



■k(a). 



Because of the wide range of extinction values and the variety 
of extinction laws appropriate for various types of galaxies (e.g. 
Calzetti 1997) we prefer to give the calibrations uncorrected 
for dust to be apphed to extinction corrected observed galaxy 
spectra. 



J. Bicker & U. Fritze - v. Alvensleben: Metallicity dependent calibrations of flux based SFR tracers 



3 



3.1. H„ 

The standard calibration for SFRs derived from H„ luminosity 
is given by K98: 

Lh„ (erg s-') = Ch„ ■ SFRh„ (Mo yr"') 

K98 derives Ch„ = 1.3 x lO"*' using population synthesis mod- 
els with solar abundances and a Salpeter IMF with an upper 
mass limit of 100 Mq. For our solar metallicity model we find 
exactly the same value for Ch„, but we also find a strong de- 
pendence of Ch„ on metallicity of more than an factor of two 
between our solar metallicity and our lowest metallicity model 
(1/50 solar) (cf. Tab.0, even though the Hq, line strength itself 
does not directly depend on the gas metalUcity. The effect is 
caused by the higher temperatures and, hence, higher ionising 
fluxes of low metallicity stars as compared to higher metal- 
licity stars. Using the standard calibration for solar abundance 
will overestimate the SFR in low metallicity galaxies by up to 
a factor > 2. 

If we choose 120 Mq as our upper mass limit, model galax- 
ies get higher Lyman continuum photon fluxes, resulting in 
stronger emission lines. Therefore, the calibration constant in- 
creases by 10 - 15%, slightly depending on metallicity, as seen 
in Tab. 1., resulting in lower SFRs at a given H,, luminosity. 

3.2. [Oil] 

The forbidden [Oil] line at 3727 A is often used as a SF indi- 
cator for galaxies at higher redshift, where the Hq, line is red- 
shifted out of the optical window. Calibrations by G-i-89 and 
Kennicutt(1992) yield 



L[Oii] (ergs ') 



C[ou] ■ SFR[oii] (Mo yr-'), 



with C[oii] = 1.5 X 10"*' and 5 x 10"*°, respectively. Both caH- 
brations are based on observed galaxy samples with SFRs pre- 
viously derived from Ho, . Differences in the galaxy samples 
are responsible for the differences in the calibration constants. 
G-i-89 used blue irregular galaxies, Kennicutt (1992) a sam- 
ple of normal, mostly spiral-type galaxies. They probably re- 
flect the lower extinctions and metallicities in irregular galax- 
ies as compared to normal spirals. K98 derived a corrected 
(IMF, H„ calibration) average for the calibration constant of 
C[oii] = 0.7 X 10**'. Our model results are close to those of 
Gh-89, which better suit to our dust-free models (cf. Tab. 
With increasing metallicity we find a rise in the calibration 
factor up to half-solar metallicity. The first step in metallic- 
ity from Z=0.0004 to 0.004 alone increases C[on] by a factor 
3. This reflects the increasing oxygen abundance. Towards so- 
lar and super-solar abundances, however, the calibration factor 
drops as a result of the lower output of ionising photons of 
high metallicity stars. In a recent work by Kewley et al. (2004) 
the metallicity dependence of the [Oil] calibration was investi- 
gated in great detail by examining the [Oil] /Ha ratio in empir- 
ical and theoretical approaches. They also find that the calibra- 
tion constant peaks at half solar metallicity, and their calibra- 
tion constants correspond fairly well to our values for Z=0.02, 
Z=0.008, and Z=0.004 with values of 1.4, 2.0, and 1.3, respec- 
tively. Unfortunately, their calibrations cannot be used for nor 



Table 2. Calibration constants for the UV luminosities vs. SFR 
relations. 



M„ 



0.0004 
0.004 
0.008 
0.02 
0.05 



100 M( 



-o 



Cl50() C28OO 
^^Q27 ei-gs-' Hz 



Mo yr- 



14.6 
12.8 
11.9 
10.6 
9.3 



13.7 
11.8 
10.8 
10.0 
8.9 



120 M, 







Ci5()0 C28OO 
27 erg s"' Hz 



xlO^ 



Mo yr- 



14.8 
13.1 
12.2 
11.0 
9.8 



13.9 
12.0 
11.0 
10.3 
9.2 



extrapolated to our lowest and highest metallicities, as their 
polynomial fit to the [Oil] /Ha vs. Z relation would then give 
unphysical negative values. 

Except for the lowest metallicity (1/50 solar), the metal- 
licity effect on the calibration of the [Oil] vs. SFR relation is 
small as compared to that of the H„ vs. SFR relation. This re- 
sults from the counteracting effects of increasing metallicity 
and decreasing number of Lyc photons. 

We reemphasize the need to accurately correct an observed 
spectrum for dust before using either H„ or [Oil] as a SFR in- 
dicator (cf. Kewley et al. 2004). 

3.3. Impact of the choice of stellar evolutionary input 
physics 

The impact of the specific choice of stellar evolutionary input 
physics can be seen from a comparison between the present re- 
sults based on recent Padova isochrones and those presented 
in Weilbacher & Fritze - v. Alvensleben (2001) using Geneva 
stellar evolution models for the three metallicities in com- 
mon, Zo, Z = 0.008, and Z = 0.004. Except for the steflar in- 
put physics, both approaches are identical, in particular, they 
both calculate the Lyman continuum photon fluxes on the ba- 
sis of Schaerer & de Koter's NLyc(Tetf) calibrations. For solar 
metallicity, H„- and [OII]-fluxes are lower by 12 and 6%, re- 
spectively, with Geneva than with Padova models, leading to 
emission line based SFR estimates higher by these percentages 
with Geneva than with Padova physics. Towards subsolar met- 
allcities, the effect changes sign, and SFR estimates from Ha 
and [Oil] are lower with Geneva than with Padova stellar evo- 
lutionary input physics by 9 and 4%, respectively, forZ= 0.008 
and by as much as 25 and 30%, respectively, for Z- 0.004. 

3.4. Luv 

The UV luminosity is also often used to derive SFRs. With the 
usable wavelength range being 1250 - 2800 A , a couple of 
calibrations are used in the literature. We here use for our cali- 
brations the wavelength ranges given by M-i-98 at 1500 A and 
2800 A which are conventionally averaged over a rectangular 
bandpass of width AA/A - 20%. 



Luv (erg s 'Hz ') 



SFRuv (Moyr-'). 
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The calibration constants obtained with our models for galaxies 
with various metallicities and different upper mass limits for 
the IMF are given in Tab.|2 In contrast to the emission line vs. 
SFR calibrations, these UV vs. SFR calibrations turn out to be 
almost insensitive to the upper mass limit of the IMF. Going 
from 100 to 120 M© only increases the calibration constants by 
less than 2 and 4% for C1500 and even less for C2800 at solar and 
low metallicities, respectively. 

Using Bruzual & Chariot (1993, BC93) models with so- 
lar metalHcity and a Salpeter IMF from 0.1 to 125 Mq, M+98 
found Cisoo = 8.0 x 10^^ and C2800 = 7.9 x 10^^ after trans- 
formation from their Ly to our L^, i.e. 27 and 23% smaller than 
our values, C1500 = H-O x 10^^ and C2800 = 10.3 x 10^^ for 
an upper mass limit of 120 M©. The difference in the upper 
mass limit has negligible effect, the diff'erence in the lower mass 
limit, however, accounts for 15% of the difference. Another dif- 
ference is that the continuous emission of the gas is included in 
our models, but not in those used by M-i-98. This causes a small 
increase in the UV luminosity for our models as compared to 
the BC93 models that M+98 use. Neglecting the continuum 
emission in our models reduces our C1500 and C2800 by another 
7%. The remaining 5 and 1% differences in €1500 and C2800 
must be due to differences in the stellar evolutionary tracks and 
model atmospheres. 

Like for Hq, vs. SFR, we find a strong metallicity depen- 
dence for the UV vs. SFR calibrations as a result of the higher 
luminosities low metallicity stars. Additionally, their higher 
Lyc photon fluxes increase the continuous emission. Taking 
all effects into account, we find that the SFR of low metallic- 
ity galaxies is overestimated by up to a factor 1 .9 if using the 
Mh-98 calibration. 



4. Conclusions 

SFR determinations for galaxies on the basis of Hq, and 
[Oil] line luminosities, and luminosities in the UV are com- 
mon techniques used to estimate the SFR in nearly all kinds of 
galaxies, from local ones to high redshifts. Even though galax- 
ies show a large scatter in metallicity, the standard calibrations 
are derived from evolutionary synthesis using solar metallic- 
ity input physics. We present here a new set of calibrations 
for the most widely used SFR indicators H„ , [Oil] , and Luv 
derived from evolutionary synthesis models and consistently 
accounting for a large range in galaxy metallicities. Using re- 
cent Padova isochrones, the spectral library from Lejeune et al. 
(1997, 1998) and recent compilations of Lyc photon rates as a 
function of stellar effective temperature (Schaerer & de Koter 
1997, Smith et al. 2002), we find good agreement with standard 
calibrations in the literature for solar metallicity galaxy models. 
At higher and lower metallicities, however, we find significant 
deviations. 

Towards lower metallicities, in particular, we find strong 
deviations from the SFR-H„ and SFR-Luv relations and pre- 
dict SFR of metal-poor galaxies derived from standard solar 
metallicity calibrations to be overestimated by factors up to 2 
due to the hotter temperatures and the higher ionising fluxes 
and UV-luminosities of low-metallicity stars. 



In case of [Oil] , on the other hand, the lower oxygen abun- 
dance and the higher output rate of Lyc photons act against 
each other and produce a maximum in Lpuj/SFR around 
Z = 1 /2 Zq. Hence, the SFRs of galaxies around half solar 
metallicity tend to be slightly underestimated while those of 
very low metallicity galaxies are overestimated by the widely 
used standard calibrations. Our results on [Oil] agree well with 
those obtained from an independent approach by Kewley et al. 
(2004). 

Our results depend on the choice of stellar evolution mod- 
els as far as hot star temperatures and lifetimes are concerned. 
Differences between Padova and Geneva stellar evolutionary 
input physics are small at half-solar to solar metallicities (< 
10%), but increase to ~ 30% towards Z- 0.004. Stellar model 
atmosphere calculations have reached a comforting agreement 
in the last years as far as the UV- and the H-ionising fluxes are 
concerned, divergences only appear around He-ionisation. 

Our results depend on the upper mass limit of the IMF with 
calibration constants for the Hq, and [Oil] vs. SFR relations in- 
creasing by 13 and 20% for low and solar metallicity, respec- 
tively, but only by less than 2 and 4% for the UV vs. SFR 
constants in the sense that SFRs assuming too low an upper 
mass limit are overestimated by these percentages. Results also 
depend on the lower mass limit, since mass normalisation re- 
quires e.g. a higher number of low mass stars be compensated 
for by a smaller number of ionising UV-bright high mass stars. 
Going from 0.15 to 0.1 M© lowers the Lyc and UV fluxes at a 
given SFR by ~ 15%. 

Spectroscopy at all redshifts picks out the most luminous 
and, hence, the most metal-rich galaxies. Hence, SFRs and SFR 
densities determined from Hq, and [Oil] fluxes are expected to 
only be affected by metallicity effects at the highest redshifts. 
SFRs and SFR densities determined on the basis of restframe 
UV-luminosities of high redshift galaxies in deep fields that 
also include the bulk of the lower luminosity and more metal- 
poor objects at each redshift, in contrast, are expected/predicted 
to be severely overestimated when using the standard calibra- 
tions. 

We therefore stress the necessity to simultaneously deter- 
mine the metallicity of a galaxy and its SFR - as can be done 
both from spectroscopy if enough emission lines are seen and 
from multi-band imaging via the comparison between observed 
and model spectral energy distributions. 
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